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Motivated by recent free-standing graphene experiments, we show how thermal fluctuations affect 
the mechanical properties of microscopically thin solid ribbons. A renormalization group analysis of 
flexural phonons reveals that elongated ribbons behave like highly anisotropic polymers, where the 
two dimensional nature of ribbons is reflected in non-trivial power law scalings of the persistence 
length and effective bending and twisting rigidities with the ribbon width. With a coarse-grained 
transfer matrix approach, we then show how thermalized ribbons respond to pulling and bending 
forces over a wide spectrum of temperatures, forces and ribbon lengths. 
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Over the last few decades the effects of thermal fluctua¬ 
tions acting on one dimensional polymers and two dimen¬ 
sional solid membranes have been studied extensively. 

It is well known that polymers remain straight only at 
short distances, while on distances larger than persis¬ 
tence length fp polymers perform a self-avoiding random 
walk mm- On the other hand, because of strong ther¬ 
mal renormalizations triggered by flexural phonons [3], 
arbitrarily large two dimensional membranes remain flat 
at low temperatures, with strongly scale-dependent en¬ 
hanced bending rigidities and reduced in-plane elastic 
constants HIS]. 

A related scaling law for the membrane structure func¬ 
tion of a solution of spectrin skeletons of red blood cells 
was checked in an ensemble-averaged sense via elegant 
X-ray and light scattering experiments [6]. However, re¬ 
cent advances in growing and isolating free-standing lay¬ 
ers of crystalline materials such as graphene, BN, WS 2 or 
M 0 S 2 Cl (not adsorbed onto a bulk substrate or stretched 
across supporting structure) hold great promise for ex¬ 
ploring how flexural modes affect the mechanical prop¬ 
erties of individual sheet polymers that are atomically 
thin. Graphene also offers the opportunity to study how 
soft flexural phonons affect the electron transport under 
various conditions [HE], and there is a prediction of a 
buckling instability in hole-doped graphene m- Experi¬ 
ments carried out in a vacuum (as opposed to membranes 
embedded in a liquid solvent) can be extended to very 
low temperatures, where the quantization of in-plane and 
flexural phonon modes becomes important muig. 

Here, inspired by recent work by Blees et al., Ha 
we consider thermal fluctuations of microscopically thin 
solid ribbons of width W and length L W. We 
show that sufficiently long ribbons behave like highly 
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anisotropic one dimensional polymers, with the two- 
dimensional nature reflected in very large renormaliza¬ 
tions of bending and twisting rigidities at the scale of 
the ribbon width W. It is natural to coarse-grain and 
construct a ribbon with L/W ^ 1 square membrane 
blocks of size WxW. Below we make this idea precise, 
by integrating out all fluctuations on scales smaller than 
the ribbons’ width. The work of Blees et al. [13] fo¬ 
cuses on the deflections and thermal fluctuations of free¬ 
standing graphene in the cantilever mode, and found a 
renormalized bending rigidity for 10/im wide ribbons at 
room temperature ^6000 times larger than its micro¬ 
scopic value at T = 0. [38] Although these ribbons were 
much shorter than the persistence length fp, which is on 
the order of meters (see below), it is possible to reach the 
semi-flexible regime (ribbon length L > ip) for narrower 
graphene nano-ribbons. With narrower free-standing rib¬ 
bons in mind, we use a coarse-grained transfer matrix 
approach to analyze the response of thermalized ribbons 
to pulling and bending for the wide spectrum of temper¬ 
atures, forces and ribbon lengths. 

To properly define the relevant quantities, we first dis¬ 
cuss thermal fluctuations of large two dimensional mem¬ 
branes under an external edge tension aij. The free en¬ 
ergy cost associated with small deformations of mem¬ 
branes around the reference flat state is El 

E = y *dxdy i [\u% + - 2hg det{Kij)] 

dr rhiaijUj, (1) 

where first two terms describe the cost of stretching, 
shearing and compressing, and the next two terms de¬ 
scribe the cost of membrane bending. The boundary in¬ 
tegral corresponds to the work done by external tension 
{fhi describes the unit normal vector in the X-Y plane 
to the membrane boundary), and summation over all in¬ 
dices i,j G {x,y} is implied. The strain tensors Uij = 
{diUj + djUi)l2 + {dif){djf)/2 and Kij = didjf, mea¬ 
sure deformations from the preferred fiat metric and zero 
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curvature respectively; we kept only the lowest orders in 
terms of the in-plane phonon deformations Ui{x^y) and 
out-of-plane deformations f{x,y) [14]. 

The effects of thermal fluctuations are reflected in cor¬ 
relation functions obtained from functional integrals jH 
E], GuiUj{r 2 -ri) = i fD[ui,f] and 

Gff{r 2 -Ti) = ^ fD[ui, /] /(r 2 )/(ri)e“'®/'=s^, where T 
is temperature, Z = fD[ui, is the partition 

function and r = (x^y). In the absence of external ten¬ 
sion {aij = 0 ), it is known that non-linear couplings of 
strain tensor Uij through the out-of-plane flexural phonon 
deformations f{x^y) produce universal power law scal¬ 
ings of correlation functions G(q) = f(d^r/A) 
in the long wavelength limit (q) rsj q ^ and 

^ where A is membrane area, y ^ 

0.82 |3l ITSHTT] and the exponents 77 ^^ + 2 ?] = 2 are con¬ 
nected via Ward identities associated with the rotational 
symmetry m- Thermal fluctuations become important 
on scales q~^ larger than thermal length [SlIlsHIHl, 

4h - Kl^/kBTY, ( 2 ) 

where Y = 4/i(/i + A)/( 2 /i + A) is Young’s modulus, 
and correlation functions above can be interpreted as 
scale dependent elastic moduli K{q)^ Kciq) ^ q~^ and 
X{q),li{q) [HIS]- Bending rigidities thus diverge 

for large membranes, while in-plane elastic constants be¬ 
come extremely small. 

In order to see the role of external tension aij 7 ^ 0, 
which will help us understand pulling forces in ribbons, 
it is convenient to integrate out the in-plane degrees of 
freedom and study ^eff = — ksT In (^fVlui] , 

the effective free energy for out-of-plane deforma¬ 
tions, HI Eeff =Jdxdy [(k/2) (VV)^ - Kadetididjf) + 
<^iji9if){djf) + {Y/8) {Plj{dif){djf))'^], where the trans- 
verse projection operator reads = 6ij — didjjY^. 
In the effective free energy description above we see 
that external tension suppresses out-of-plane fluctua¬ 
tions in /, which have long range anharmonic interac¬ 
tions between transverse tilt deformations of the mem¬ 
brane normals. The effects of the anharmonic term at 
a given scale t = 27r/g* can be obtained by integrating 
out all degrees of freedom on smaller scales. Formally 
this is done by splitting all fields ^(r) G /(r)} 

into slow modes 5 '<(r) = 

modes ^>(r) = X]|q|>g* which are then inte¬ 
grated out as E{JP) = —/c^Tln />] . 

The functional integrals following from standard per¬ 
turbative renormalization group calculations na na 
uni lead to a free energy with the same form as 
in Eq. 0 except that renormalized elastic constants 
\R{t),iiR{P),nR{P),t<iGR{P) become scale dependent, 
while the external tension a^j remains intact. 

For a small isotropic external tension aij = a5ij , or for 
a small uniaxial tension in the x-direction aij = aSix6jx^ 


the tension becomes relevant on scales larger than [20] 


f K 

\<J 


(3) 


where exponent y ^ 0.82 and thermal length scale ith [see 
Eq. 0] have been defined above for membranes without 
external tension. As shown in the supplemental mate¬ 
rials, external tension then produces the renormalized 
elastic constants 


i^r{^) i^rg{^) 

R ’ RG 

^r{^) 

A ’ y 


1 , ^ < 4h 

W4h)^ 4h<^<4 , 
(4/4h)^ 4 < ^ 

1 , ^ < 4h 

4h<^<4 .(4) 


The out of plane correlation function then becomes 
(q) = ^ [KR{2'K/q) q'^ + aijqiqj] . For isotropic ex- 
ternal tensions this result agrees with Roldan et al [20] , 
but the results for uniaxial external tension appear to 
be new. With a uniaxial tension, the long wave length 
/(q) fluctuations behave like the layer displacements of 
a defect-free two dimensional smectic liquid crystal [ 21 ], 
with fluctuations along the direction x of the pulling 
force having a reduced amplitude Gjj (|q| < ^ 

^ [crg2 + Kgj^4(4/4h)’'] • 

For sufficiently large external tension a > ksTY/R = 
cr*, [ 20 ] which corresponds to 4 h ^ thermal fluc¬ 
tuations become irrelevant and the renormalized elas¬ 
tic constants are approximately equal to the microscopic 
ones. Remarkably, for graphene membranes with r = 
l.leV [22] and Y = 340N/m [23], the thermal length 
at room temperature is of order the lattice constant, 
ith ^ 1^! [SI US] Therefore thermal fluctuations are im¬ 
portant for all experimental situations in this case pro¬ 
vided only that the external membrane tension is smaller 
than cr* ^ lON/m. 

We now study ribbons of width W and length L that 
lie on average in the X-Y plane with long axis in x in 
direction and with a pulling force F = Waxx on the 
ribbon end. Once we integrate out all degrees of freedom 
on scales smaller than VF, the resulting strain tensors 
Uij and Kij depend only on the x coordinate, and the 
renormalized elastic constants in Eq. 0 are evaluated 
at = W. This results in an effectively one dimensional 
free energy model for the ribbon 

W 

-£'= / da; Y [Xru% + 2fiRujj + krK^ - 2k.gr det{Kij)] 

-Fux{L). (5) 


If we then continue integrating out degrees of freedom 
on scales larger than W in this effective one dimen¬ 
sional problem (see supplemental materials), the renor¬ 
malized bending moduli rr^rgr and the renormalized 
shear modulus jiR remain constant. However, the in¬ 
plane elastic modulus 2yR + Xr, which is related to the 
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FIG. 1: (Color online) (a) Ribbon configurations with unde¬ 
formed length L > W can be described with orientations of 
material frame { 01 , 02 , 63 } attached to the ribbon relative to 
the fixed laboratory frame {x,y,z}. (b-c) Initial ribbon ori¬ 
entations for studying the response (z) to external bending 
and pulling forces F. 

deformations Ux{x,y) averaged over the ^-direction, be¬ 
comes smaller and smaller, a sign that the ribbon does 
not remain straight. In this regime the free energy de¬ 
scription of small deformations around the flat state in 
Eq. <§ breaks down for small F. 

For long ribbons L > W we exploit a complementary 
description that allows for large deformations in three 
dimensional space, provided that local strains remain 
small. We attach a material frame {ei(s), 62 ( 5 ), 63 ( 5 )} 
to the ribbon (see Fig. [^, where s G [0,1/] is the co¬ 
ordinate along the center of the ribbon backbone. The 
orientation of material frame relative to the fixed labo¬ 
ratory frame {x,y,z} can be described with Euler an¬ 
gles S{s) = {0(s), 6 ^( 5 ),'0(s)} [24]. The laboratory 
frame can be rotated to the local material frame with 
three successive three-dimensional rotations R(©) = 
Rz{—i^)Ry{—0)Rz{—(t)) [25l|26], such that { 01 , 62 , 63 } = 
{R(©)x,R(©)y,R(©)z}. 

Ribbon bending and twisting deformations are then de¬ 
scribed by the rate of rotation of the material frame along 
the ribbon backbone ^ = ^R“^6q, = ft x [25l \T7\ 
with ft{s) = Here and ^^le radii 

of curvatures for bending of ribbon around axes 61 (s) and 
62 ( 5 ), and 27rO^^(s) describes the pitch for ribbon twist¬ 
ing. The free energy cost of a ribbon deformation is then 
written as [251427] 

E= /ds 1 + ^2^1 + cnl] - F ■ r(L), (6) 

Jo 2 

where Hi, H 2 are bending rigidities, C is torsional rigidity 
and F is the applied force on a ribbon end at which 
can be a bending or pulling force (see Fig. [H. From 
comparison with the effective one dimensional model for 
ribbon in Eq. we find 

= Wi^r{W), C = 2Wi^gr(W), (7) 

where renormalized bending rigidities are defined in 
Eq. (|4|. The second bending rigidity for bending around 
axis 62(5), involves ribbon stretching and is much larger; 
in fact, H 2 ’s bare value exceeds Hi and (7 by a fac¬ 
tor of order YW^ j the Foppl-von Karman number. 


where we expect kq ^ R for graphene both microscop¬ 
ically El and when thermal renormalizations are ac¬ 
counted for. This quantity can be estimated from clas¬ 
sical zero temperature solid mechanics m as H 2 ^ 
W^Yji{W). Here, the renormalized Young’s modulus 
Yr = 4/ii?(/ii? + \R)/{2fiR + \r), scales in the same way 
as other in-plane elastic constants \r and fiR in Eq. 

For ribbons whose width is much larger than it’s thick¬ 
ness we thus find H 2 ^ Hi,(7 and we can set Vt 2 ~ 0 . 
In the free energy cost for ribbon deformations [Eq. Q], 
we neglected the stretching/compressing of ribbon back¬ 
bone, as is appropriate when the pulling force resisting 
entropic contraction is not too large [28]. The effective 
one dimensional free energy model presented above cor¬ 
responds to the highly asymmetric Id polymer [251427] . 
with anomalous kF-dependent elastic parameters. [39] 

The free energy described by Eq. is complicated 
because successive rotations do not commute! However, 
the physics can be understood by mapping the statisti¬ 
cal mechanical problem to the corresponding quantum 
mechanical problem [29] as described below. 

The response {z) of the ribbon to external force F in 
the z direction can be evaluated from the relation {z) = 
ksT^dhiZ/dF)^ where the partition function reads Z = 
/P[©(s)] e ^ Note that we can study both pulling 

and bending forces, where the only difference is the con¬ 
ditions on the Euler angles, i.e. in the initial orientation 
of ribbon (see Fig.Q. If we clamp the ribbon at the ori¬ 
gin (5 = 0) and apply force on the ribbon end {s = L), 
then for pulling the initial condition is ©^ = {0,0,0}. 
To treat bending, we consider a ribbon initially aligned 
with the x-axis and take ©^ = {7r/2,7r/2, 0}. To eval¬ 
uate the partition function Z, it is convenient to define 
the unnormalized probability distribution p(©, s) of Eu¬ 
ler angles © at a contour length s along the ribbon mid- 
line as p{®f,Sf) = , where 

the path integral above is restricted to s G [0,s/] and 
the partition function is given hy Z = fdSp(S, L) with 
the Euler-angle measure fdS = d(j) sin Odd 
The evolution of this probability distribution along the 
ribbon backbone is described with differential equation 

~ O 9 [251 l3Qj where the Hamiltonian 

operator is .ff = ^ (^ + ^ + #) - Here 

the {Ja} are angular momentum operators around axes 
6 q,, which can be expressed in terms of derivatives with 
respect to Euler angles PlIHUj. The evolution of p(©, s) 
with s maps the physics of thermalized ribbons onto the 
Schrodinger equation of the asymmetric rotating top m 
in an external gravitational field, where the ribbon back¬ 
bone coordinate s plays a role of imaginary time and 
the bending and twisting rigidities Hi, H 2 and C cor¬ 
respond to moments of inertia. The evolution of the 
material frame orientation distribution can be evaluated 
by expanding the initial condition in eigen-distributions, 
p(&,0) = 6(& - ©i) = T,aBaPa{®), where Hpa{&) = 
^apa{^)- In this decomposition the partition function 
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FIG. 2: (Color online) Pulling and bending deflections {zjL) 
of ribbons with bending rigidities A^jAx oo and twisting 
rigidity CjA\ — 1 in response to a fixed small external force 
FAll{ksTY — 0.01. The slope of +2 for bending when L 
Ip agrees with expectations for stiff cantilevers with, however, 
a bending rigidity greatly enhanced by a factor (VF/Gh)^ ^ 
1. The responses to pulling and bending forces agree when 

T >4. 

becomes Z = f dSpa(S) and the response 

(z) to an external force can be evaluated from Z as de¬ 
scribed above. 

To treat ribbons in both the semifiexible and highly 
crumpled regimes, we must now find all eigenvalues 
and eigen-distributions Pa(^)- From quantum mechan¬ 
ics we know that this is done efficiently in the basis of 
Wigner D functions [24], which have well de¬ 

fined quantum numbers j^k^m for the total angular mo¬ 
mentum + JI + J^ 3 , the angular momentum 

around the ribbon tangent J 3 and for the angular mo¬ 
mentum around the laboratory axis Jz- For details see 
Refs. [25l |30| and the supplemental materials. 

With the help of this machinery we first studied the 
response of ribbons of various lengths to small exter¬ 
nal pulling and bending forces at fixed temperature (see 
Fig.§. Here, since C and Ai have a similar order of mag¬ 
nitude, we take C = Hi, for simplicity. Similar to single 
molecule polymer physics, [311132] we find two regimes. 
For ribbons much shorter than a persistence length [26] 

/ 2 2WKRiW) 

^ kBT{Ap+Ap)^ keT ■ 

ribbons behave like stiff “classical rods” El, where for 
pulling {z) ^ L and for the bending (cantilever) mode 
{z) = FL^/3Ai. Note that A 2 ^ is negligible and that 
thermal fluctuations on scales less than W lead to a 
renormalized bending rigidity Hi [see Eq. 0 ] , orders 
of magnitude larger than for rod-like polymers at room 
temperature, as found by the Cornell experiments [T3] . 
For ribbons much longer than the persistence length 
{L ^ ip)^ pulling and bending become equivalent. In 
this semi-flexible regime ribbon forgets its initial orien¬ 
tation after a persistence length, and for small pulling 
forces the response to either bending or pulling is {z/L) = 



FIG. 3: (Golor online) Response of ribbons (neglecting quan¬ 
tum fluctuations) to a small bending force at various temper¬ 
atures for fixed VF, L, F, k and Y. Three regimes appear 
for the parameter choices, FI? jZW n — 0 . 01 , jn — 10 ^, 

L/FF = 10^ 


2Fipj(?ikBT^ [2]- Eventually, at much larger ribbon 
lengths than those considered here, ribbon self-avoidance 
will become important. [28] 

To highlight the difference between conventional poly¬ 
mers and thermalized ribbons with FF ^ ^th, con¬ 
sider the response of ribbons to a small bending force, 
FL^/Wk < 1. Eig. I shows results for a wide vari¬ 
ety of temperatures, obtained by inserting temperature 
dependences hidden in Hi and ip. We find three dis¬ 
tinct regimes: At small temperatures, where FF ^ ^ 

^/ksTY^ thermal fluctuations are negligible and rib¬ 
bon behaves like a classic cantilever with bare elastic 
parameters, ( 2 :) = FL^I{3kW). As the temperature in¬ 
creases, the thermal length scale drops and eventually 
becomes smaller than the ribbon width (^th ^ kF)- In 
this regime the renormalized bending rigidity is increased 
due to thermal fluctuations and the cantilever deflection 
is smaller {z) = FL^i^y^/{3hzW^~^^) ^ As tem¬ 

perature increases even further, eventually the persis¬ 
tence length ip becomes smaller than the ribbon length 
L. As noted above in this semi-flexible regime the de¬ 
flection now becomes {z) = 4 :kFLW^~^^/{ 3{kBT)‘^i^^) ^ 
drops even faster with temperature, as 
the ribbon transforms from a cantilever into a random 
coil. Note that with rising temperatures the cutoff length 
scale i(j associated with ribbon tension [see Eq. 0 ] also 
increases, but never becomes relevant. 

However, ribbons with large pulling forces nevertheless 
show a non-trivial response due to the cutoff i^ • Eor large 
pulling forces, Fip ^ kBT^ we also need to include the 
stretching of the ribbon backbone, with the result similar 
to Ref. [25] (I;) « 1 -h (see also supplemen¬ 

tal materials), where Yid = WYb{W) is the effective one 
dimensional ribbon Young’s modulus. The middle term 
describes stretching of the ribbon backbone, and the final 
correction corresponds to the entropic contribution from 
ribbon fluctuations. As F = axxW increases the cutoff 
length scale ia- [Eq. Q] drops and we find two crossovers. 
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FK./{kBTYW) 


FIG. 4: (Color online) Contributions of backbone stretching 
(blue line) and entropic elasticity (dashed black line) describ¬ 
ing the response to large ribbon pulling forces. We chose 
parameters ksTj k, = 1/40 (suitable for graphene at room 
temperature), and Wjith = 10^ (1 /im width ribbon at room 
temperature). 


first when this length scale crosses the ribbon width W 
and finally when it drops below the thermal length scale 
Ah (see Fig. i- Especially interesting is the intermediate 
force regime with £th ^ W, where we find that the 

ribbon backbone stretches as F/Yid ^ ^ which 

generalizes to ribbons the result of Ref. [16] for the non¬ 
linear stretching of two dimensional membranes under a 
uniform tension. See supplemental materials for further 
discussion. 

For graphene ribbons, where ^th ^ lA at room tem¬ 
perature, the experiments of Blees et al. Ha on ribbons 
of width W = 10/im, confirmed a renormalized bend¬ 
ing rigidity nR{W)/n ^ 6000, consistent with Eq. Q. 
The corresponding persistence length is of order of me¬ 
ters. Thus ip ^ L for graphene ribbons of lengths 
L ^ 10 — lOO/rm, which should behave like conventional 
cantilevers with, however, a strongly renormalized L in¬ 
dependent bending rigidity. Probing the semi-flexible 
regime requires narrower ribbons of order 10 nanometers 
width, so that the persistence length should be in the 
experimentally accessible regime of 10-100 micrometers. 
Although the value of the critical pulling tension, beyond 
which thermal fluctuations become irrelevant, is F/W = 


cr* ^ lON/m for graphene, one could observe interesting 
behavior for smaller tensions where it^ < i^ < W. Ad¬ 
ditional novel behavior can arise for free-standing sheets 
at sufficiently high temperatures even when L ^ W. To 
see this, consider the correlation function of the mem- 
brane normals n{x,y) = {-d^f, -dyf, 1)/^/! + |V/p 
that defines the flat phase [3]. There is a power law 
approach to long range order, (n(ra) • = 1 ~ 

+ln(Ah/ao)] + where C is 

a positive constant of order unity and ao is microscopic 
cutoff, of order the graphene lattice spacing (see supple¬ 
mental materials which includes the effect of an isotropic 
external stress). The second term represents the reduc¬ 
tion in the long range order due to thermal fluctuations. 
When this term becomes the same size as the first (i.e. 
for ksT > 271KT])^ the low temperature flat phase should 
transform into a entropically dominated crumpled ball, 
with a size limited by self-avoidance, provided monolayer 
sheets such as graphene maintain their integrity [18] . The 
transition temperature to isotropic crumpling could be 
lowered by creating a graphene sheet with a periodic ar¬ 
ray of holes or cuts. (Although cuts could be deployed 
with equal numbers at 120 degree angles, an array of par¬ 
allel cuts could lead to a system that is crumpled in one 
direction, but tube-like in another, a situation studied 
theoretically in Ref. [33].) While we have some under¬ 
standing of force-free conformations [34], little is known 
about the mechanical properties of free-standing mem¬ 
branes at or above this crumpling transition. There is 
evidence from computer simulations of a high tempera¬ 
ture compact phase, where attractive van der Waals in¬ 
teractions are balanced by self-avoidance [35]. We hope 
this paper will stimulate further investigations on these 
problems in the spirit of single-molecule experiments on 
linear polymers [snisi]. 
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Supplementary Information 


In this supplemental material we provide detailed calculations that were omitted in the main text for clarity. Section 
I describes the renormalization group treatment of membranes under tension. In Sec. II we present how pulling and 
bending of ribbons can be mapped to a time evolution of rotating top in quantum mechanics, and how tools from 
quantum mechanics can be used to calculate the force-extension curves for ribbons. 


I. RENORMALIZATION GROUP TREATMENT OF MEMBRANES UNDER TENSION 


Our goal is to analyze properties of fluctuating membranes under external tension cjij with the renormalization 
group approach. The free energy cost of membrane deformations under tension is 


E = 


J i [Xu% + 2jj.ufj + - 2k,g det{Kij)] - ^ds 


ElidijVjj , 


(SI) 


where Uij = {diUj + djUi + difdjf)/2 is the nonlinear strain tensor, Kij = didjf is the bending strain tensor, the 
Ui are in-plane deformations, / is the out-of-plane deformation, and mi describes a normal vector to the membrane 
boundary. Using the divergence theorem we can convert the boundary work term to the area integral, such that the 
free energy becomes 


E = J Q [Xu% + 2ijujj + kK% - 2kg det{Kij)] - , (S2) 

where = {diUj + djUi)l2 is the linear part of the strain tensor. Since the in-plane deformations Ui only ap¬ 
pears quadratically in (Eq. [S^ , we can integrate them out to derive the effective free energy for the out-of-plane 
deformations, [4] 


E 

~A 


E \ /(q)/(-q) 

q 

+ E ^ [9ii-Fiyq)92i] [g3i-F)yq)?4i] /(qi)/(q2)/(q3)/(q4), 

qi+q2=q7^0 

q3+q4=-q7^o 


(S3) 


where the Young’s modulus is Y = -f A)/(2/i -]- A), the projection operator PIM) = % — qiqjlq^^ A is the 

membrane area and the Fourier modes are /(q) = / {d?Y/f {y). From the expression above we can clearly see 
that positive components of the membrane tension dij constrain the out-of-plane fluctuations /. 

To implement a momentum shell renormalization group, we first integrate out all Fourier modes in a thin momentum 
shell A/h < q < where A is microscopic cutoff and h = iA = with d ^ 1. Next we rescale lengths and fields mm 

X = 6x', 

/(x) = 6^//'(x'). (S4) 


We find it convenient to work directly with a D = 2 dimensional membrane embedded in d = 3 space, rather than 
introducing an expansion in e = 4 — D. m Finally, we define new elastic constants Y', and external tension crF, 
such that the free energy functional in Eq. (S3) retains the same form after the first two renormalization group steps. 
It is common to introduce (3 functions [19], which define the flow of elastic constants 


P. 

Py 

Pij 


dn' 

dlnb 

dY' 

d\nb 


— 2(C/ ~ 1)^ + 

= 2(2C/ - 1)Y + Yy, 


9 In 6 




(S5) 


Above we introduced Z functions, which result from the integrals of modes over the momentum shell. To one loop 







FIG. SI: One loop corrections to the renormalization of (a) Hi and (b) Y. Solid lines represent propagators for the out-of-plane 
displacements /(q) and dashed lines represent the momentum carried by the vertex Y. 


order (see Fig. SI), the Z functions read 


d 


F (i-(q-p)") (/(p)/(-p)) b 


f <p<A 


= - 77 ^ ( ;b4; Y (l-(q-P)^)^P^(/(P)/(-P))^ 


d\nh \ 2kBT 


(S6) 


^<P<A 


where A is the microscopic momentum cutoff and A is the undeformed membrane area. Note that the only change 
in the stress tensor (Jij to this order arises from the rescaling factor Q. Upon assuming that the initial membrane 
tension (Jij is small, such that cjij <C nh? ^ then {f{p)f{—p)) ~ ksT/ (Anp^) in equations above and the P functions 
in one loop approximation become 


Pij 


2 (C/ -1)^ + 


SYksT 


2(2C/ - 1)U 


SY^ksT 




It is convenient to chose (f such that 13^ = 0, which results in (f = 1 — and 


Py 

Pij 


2Y - 


2 1 


15Y‘^kBT 
‘32TTti?A? 
SYksT \ 
327r/^^A^ / 




(S7) 


(S8) 


By repeating the renormalization group procedure, we integrate out modes at the smallest length scale and evolve 
the Young’s modulus Y and external tension aij. Initially, they both grow rapidly 

Y(i) ^ Y X (£Af, 

Gij(f) « X (M)^, (S9) 


where we integrated out all modes on scales smaller than £. Once we integrate out all modes up to the scale 
^th ^ KjyJkBTY Young’s modulus reaches a fixed point 




IbkeT 


- Yx (4hA)^ 


(SIO) 


At the fixed point we introduce the exponent r^, such that Q = 1 — r]/2. Note that C/ ~ 1 initially, before we reach 
the fixed point. In the one loop approximation we find r] = 4/5, which approximates the value of r] ^ 0.82 obtained 
by the self-consistent screening approximation m and T] ~ 0.85 obtained by the non-perturbative renormalization 
group calculations m- This result differs from a formal one loop e = A — D expansion, which results in 77 = 
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12e/25, [15] because we have performed the one loop calculations directly for D = 2 dimensional membranes, rather 
than calculating them for small e, i.e. for D ^ 4 dimensional membranes. 

By continuing with the renormalization group procedure and integrating out modes beyond the scale we find 
that the initially small membrane tension now grows as 

> 4h) = cTy X X (4hA)^- (Sll) 


Eventually, the membrane tension becomes large enough that it becomes important. This happens at the scale 

l/(2-r7) 




(S12) 


when (Tij{i(j 


. At this stage, we have to take into account the membrane tension, when evaluating the Z 
functions in Eq. (S6). In next subsections, we describe what happens for membranes under various external tension 
conditions. We first discuss membranes with W ^ L, and then move on to discuss ribbons with L ^ W. 


A. Membranes under uniform tension 


We first consider membranes under uniform tension (Jij = (j5ij. After integrating modes on scales smaller than 
the membrane tension becomes relevant and beyond this point we can approximate {f{p)f{—p)) ~ kBT/{Aap‘^) in 
Eqs. ( [S^ . With this change, the [3 functions become 


Pij 


2(C/ -1)^ + 


SYksT 

IGttct 


2(2C/ - 1)T 


SY^ksTA^ 

327rcr^ 




(S13) 


It is now convenient to set = 0 so that the uniform tension remains unchanged. We then find that both the bending 
rigidity hz and the Young’s modulus Y flow to 0 at large length scales, 

k{£ > ia) ~ K X 

Y(i > 4) ~ y* X (S14) 


In this regime the external tension dominates and thermal fluctuations are unimportant. After rescaling lengths and 
fields back to the initial units we find that the height correlation function for the out-of-plane flexural phonons is 


(/(q)/(-q)) 


kBTYAnq^), 9 > 

ksTYAaq^), £P > q 


keT 

A{aq^ + KR{q)q'^) 


(S15) 


It is useful to calculate how the membrane expands due to the external tension in the presence of thermal fluctuations. 
We find that the fractional area change is 


SA 

~A 

SA 

~A 


(4> = 


keT 


(/i + A) 


(/(q)/(-q)), 


—h In(^thA) 
A 


ksT 

Attk, 


-1 


ksTYj 


a 

{p + A) 


(S16) 


The first term above corresponds to the shrinking of membrane due to thermal fluctuations, and reflects a negative 
coefficient of thermal expansion when cr = 0, 


I dA 

A^ 


ks 

Alin 


-X + ln(4hA) 

f] 2 


(S17) 


The second term describes the nonlinear stretching for small uniform tension a in the presence of thermal fluctua¬ 
tions [T6|. The last term is the conventional linear response result, describing membrane stretching in the absence of 
thermal fluctuations. This last term only dominates for large tensions, a > ksTYlhc^ which corresponds to ^th ^ 
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Finally, we present the correlation function of the membrane normals h(x,^) = (—9^,/, —dyf, 1 )/^/! + |V/p that 
defines the flat phase [3] . When deformations are small the correlation function of the membrane normals is approxi¬ 
mately 


(n(ra) ■ n(r6)) 


1-E- 


I _ giq-Cra-rb) 


(l/(q)P). 


(S18) 


For small tension a < ksTY /k this correlation function evaluates to 
(n(ra) • n(rb)) 


1 - [V + ln(4hA)] + 77 ^( 7 ?-^ - 2-1) ' ' 


{271^) 

ksT ) 


{27171) 

4h < |i*a - r^l < 4 
y 4 < \ra - r^l 


ksTYj 


(S19) 


where C = ^ ~ 0 . 2 , D is another constant of order unity and Jo{x) is the Bessel function of the first 

kind. The second term in the equation above represents the reduction in the long range order between normals due to 
thermal fluctuations and the third term shows how this long range order is restored with external tension. For large 
tension a > ksTY /where the effects of thermal fluctuations are suppressed, we find 


(n(r„) ■ n(rb)) 1 


ksT 

--—- In IH- 

{AtTK) L (7 



(S 20 ) 


where Ko{x) is the modified Bessel function of the second kind, which asymptotically scales as Ko{x) x y^Tr/{2x)e~^. 

Note that the nonlinear dependence of the membrane extension (u^i) on the external tension a can be obtained 
simply from the scaling arguments. Since the external tension a is a conjugate variable to djUi^ their rescalings are 
connected. Once we rescale lengths as x = hx' and in-plane deformations ui = then the external tension rescales 

as cr = b^^a' with 4 = 1 ~ ^ ~ Cu? where D = 2 is the membrane dimensionality. We also know that the Ward 
identities associated with rotational symmetry connect rescaling of the in-plane and out-of-plane deformations such 
that 4 = 24 — 1 and therefore 4 = ~‘^Cf [IS]- mentioned above we can extract exponent 77 from Q = 1 — 77 / 2 , 
which leads to 4 = 1 — ^ and 4 = -2-1- 77 . Now we have all necessary ingredients to calculate the scaling of membrane 
extension as 

. (S21) 

Since the rescaling factor b is arbitrary, we can pick b = to find 

(^6u^^{a)) = (^6u^^{l)) = const, x (S22) 

Thus we found the same nonlinear scaling between the membrane stretching and the uniform tension as in Eq. 
which holds for small uniform tension. 


B. Membranes under uniaxial tension 


In this section we consider membranes under uniaxial tension axx > 0, while ayy = a^y = 0. Upon again integrating 
out modes on scales smaller than the membrane tension becomes important and beyond this point we have to 


take {f{p)f{—p)) ~ kBT/[A{crxxPx + ^Py)] Eqs. (S 6 ). Although we can ignore a term 7v{Px + ‘^PxPy) compared to 


^xxPx^ we have to keep the term with npy. Once we integrate out modes from a thin shell A/b < p < we find that 
the quadratic term in the free energy becomes 


1 
2 


/(q)/(-q) i (^xxqI + + 


YkeTlnb 

4:71 


JSL + 1 (. 


■Mt + QqIiI + 9y) 


(S23) 


All new generated terms that involve Qx are negligible compared to the ctxxQx- Therefore we can keep only the last 
term with to calculate the (3^ function that renormalizes the bending rigidity, 




2(C/-1)« + 


YkeT 

4707 


(S24) 
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For the quartic term with momentum-dependent Young’s modulus Y we find that after the momentum shell integration 
there are again anisotropic contributions in terms of and qy. Significantly, all renormalizations of Y are negative 
and (3 functions now take the form 


p. 

Py 

Pij 


2 (C/ - + 


YkeT 

A.'ko ’ 


2(2C/-l)r-|Zy|, 




(S25) 


As for the uniform tension case, we choose C/ = 0 to fix the uniaxial tension, and find that both the bending rigidity 
K and the Young’s modulus Y again flow to 0 as 

k{£ > £a) ~ K X 

Y(£ > 4) ~ Y* X (44)“^- (S26) 


After rescaling lengths and fields back to the initial units we find that the correlation function for the out-of-plane 
deformations becomes highly anisotropic 


(/(q)/(-q)) = 


kBTj{AKq^), 


q > 4 


nh 

C > 9 > 


1-1 = 


kBT 


^xxQx + KR{q)q‘^)' 


(S27) 


kBT/{A[axxql +K4(4/4h)'']), 4^ > 9 
We can now use this result to calculate the membrane strains associated with uniaxial stretching 

^xx 1 




2 ^ 
q 


«> 

«> 


ksT r _i , . x-| ksT r 

[t] ^ + ln(4hA)] + - ^ 


^^xa 


Sttk . 


f] 


— 1 + a/2 — sinh ^(1) 


tzax 


\ {ksTY) 


^77/(2-77) ^ 


XX 


ksT 


Sttk, 


“ ^ (/(q)/(-q)) > 

ln( 4 hA)] 


ksT 

Sttk. 


-1 


+ 1 — \/2 — sinh ^(1) 


/ KCTxx \ 

\k^J 


77/(2-77) 


~Y~ 


(S28) 


where z/ = A/(2/i+A) is the two-dimensional Poisson ratio. Again, the first terms in the second and fourth lines describe 
membrane shrinkage due to thermal fluctuations, the second terms correspond to nonlinear membrane stretching in 
the presence of thermal fluctuations, and the last terms correspond to the zero temperature response, which becomes 
relevant for Gxx ^ ksTY/n. The power law scalings above are accurate, but the numerical prefactors are approximate. 
In order to calculate numerical prefactors exactly, we would need to know how the correlation function in Eq. (S27) 
behaves in transition regions. In principle, the renormalized Poisson’s ratio is calculated as 


= 


(Kv) 


(S29) 


where {3u^j) describes the r elativ e change, when the uniaxial tension is increased from zero to cFxx- Because our 
numerical prefactors in Eqs. (S28) are just approximate we cannot determine the precise value of the renormalized 
Poisson’s ratio vr in the regime dominated by thermal fluctuations, but we know that the vr transitions to the zero 
temperature value v for large pulling tension, i.e. cFxx ^ ksTYj r. 


C. Pulling of ribbons 

Einally, we comment on pulling on large aspect ratio ribbons of length L and width W L. After integrating out 
all degrees of freedom on scales smaller than the width IF, the resulting strain tensors Uij and Kij depend only on 
the X coordinate and the renormalized elastic constants are evaluated at q = 27rjW. This results in an effectively one 
dimensional free energy model for the ribbon 

^ = J dxW [>iRUii + + R-rKu - 2KGRdet{Kij)] - ■ (S30) 
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P 



FIG. S2: One loop corrections to the renormalization of 2/i + A. Here, the solid and wiggly lines represent propagators for the 
out-of-plane displacement f(q) and for the in-plane displacement Ux{q), respectively. 


It is convenient to rewrite the effective free energy above in terms of one-dimensional Fourier variables s{q) = 
J{dx/s{x) and to separate out the uniform strain [4] The resulting free energy reads 

^ ^ Afl {u%)^ + (u°F + I T] + ^R)(f\ux{q)\^ + I^Rq^luyiq)^] - cr^xU^xx 


WL 


2 L 


9 T! q^\f{q)? + ^ X] + >^R)qiq2{qi + q2)ux{qi)f{q2)f{-qi - 92) 


< 71,<72 


(2/Xfl +Afl) ^ qiq2q3{qi + q2 + q3)f{qi)f{q2)f{q3)f{-qi-q2-q3)- 


(S31) 


Because the in-plane deformations Uy{q) decouple the shear modulus jur does not get further renormalized. Similarly, 
we find that the bending rigidity kr does not get further renormalized. To see this, we integrate out the in-plane 
modes {u^j^Ui{q)} to derive the effective free energy 




(S32) 


However, the in-plane modulus 2iir + \r associated with the in-plane deformation Ux{q) suffers significant renormal¬ 
izations. This can be shown with the momentum shell renormalization group by integrating out all Fourier modes in 
a thin momentum shell A/b < q < A and rescaling lengths and fields as 

X = hx\ 

Ui{x) = 

f{x) = b‘^ff{x'), 

^XX 


g' 


(S33) 

Note that the momentum cutoff is now A = 27: jW^ because we already integrated out all degrees of freedom on scales 
smaller than W. As in previous sections we define P functions that dictate the flow of elastic constants 


/?« 

/^2/x-|-A 

Pa 


dn' 

9 In 6 
dji' 

9 In 6 

9(2/i + A)^ 


9 In 6 


= 2(C/ - 1)«, 

‘^Cup’) 

= 2Cu(2/i + A) — 


dc^xx 
9 In 6 


— Ca^x 


(S34) 


The Ward identities associated with rotational symmetry connect rescaling of the in-plane and out-of-plane deforma¬ 
tions such that Cu = — 1 [16] and ("cr = — 1 — because Gxx and dxUx are conjugate variables. The integrals of 

modes over the momentum shell now only affect the in-plane modulus 2/i + A and to one loop order (see Fig. S2) we 
And 


d 


^ 2 /i+A 


9 In 6 


E pHf{p)f{-p)? 


f<P<A 


(S35) 
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Upon assuming that the external tension is small, i.e. axx ^ and choosing ("j = 1 to fix the bending rigidity 

the flow of elastic constants is described by the ;d-functions, 

= 2(2„ + A)-<'|±|^, (S36«) 

Pa- = (S36b) 


Note that in the equations above the width W also gets rescaled according to lU ^ W/{£A). If the in-plane modulus 
2 /i + A was small, then we would expect it to grow as 

+ X{£) ~ {2fiR + Xr) X {£Af. (S37) 

This modulus would keep growing until we integrate out all degrees of freedom up to the scale 


t 


^ksTipifiR + \r) 




1/3 


(S38) 


where the second term in the P 2 ^+x function in Eq. (S36a) becomes relevant. However, for small tension Gxx ^ 
we find ^ lU, because the elas tic mo duli above have already suffered large renormalizations out to the scale W. 
Therefore the second term in Eq. ( S36a| ) has to be taken into account immediately and the in-plane modulus flows as 


2/i(^) + A(^) - (2/i^ + A^) X {W/£) 


(S39) 


ioi £ ^ W. This modulus keeps dropping until the external tension becomes relevant at scale £a- ^ As 

in previous subsections the external tension introduces a cut-off length scale for the renormalization of the elastic 
modulus. By rescaling lengths and fields back to the original units we find the in-plane correlation function of 
displacements along the ribbon axis. 


{ux{q)ux{-q)) 


ksT f l/(q^W^), 

LW(2fiR + Xr) [ 4i(g^W^), 4 « 


(S40) 


Note that for small external tension the renormalization produce s larg e in-plane fluctuations suggesting that the 
description for the effective one dimensional free energy in Eq. (S30), which assumes small deformations about an 
approximately flat ribbon geometry, must eventually break down. In the next section we discuss how to treat ribbons 
with large deformations. 


II. FORCE-EXTENSION CURVE OF RIBBONS DUE TO THERMAL FLUCTUATIONS 

Consider a long thin ribbon of length L, thickness h (atomically thin for graphene!) and width W in which we 
embed a position-dependent orthonormal triad frame {ei(s), 62 ( 5 ), 63 ( 5 )}. Here s G [0,1/] is an arclength coordinate 
along the ribbon midline, 63 is a unit tangent vector along this backbone, and ei and 62 are unit normal vectors to 
the backbone as sketched below. 


e 2 



One way to express the rotation matrix R, which rotates the fixed laboratory frame {x, y, z} to the ribbon frame is to 
use Euler angles © = {0, ^,'0}, [24] via the decomposition R(S) = Rz{—'ip)Ry{—0)Rz{—(l)), such that { 01 , 62 , 63 } = 
{R(0)x,R(0)y,i?(0)z}. Here 


Ry{OL) — 


— sma 

0 


and Rz ((a) = 


cos Of 0 
0 I 
sin a 0 


cos a 


cos a sin (a 0 
— sin a cos a 0 
0 0 I 


(S41) 
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matrices correspond to rotations around the fixed laboratory axes y and z. 

Because of the rotational and translational invariance, the free energy cost of ribbon deformations only depends on 
derivatives of the attached frame, which can be expressed as the rate of rotation ft{s) of the ribbon frame [25l [26] 


dej 

ds 


(ii? 1 


ei = nx Bi, 


i = l,2,3. 


(S42) 


Here, deviation from flatness is measured by the components of ft = fliOi^ where fli(s) and fl 2 {s) are the ribbon 
bending curvatures around axes ei(s) and 62 ( 5 ), and ft 3 {s) is a twisting strain of the ribbon around the 63 ( 5 ) axis. 
Alternatively we can view fti{s) as the rates of rotation of the ribbon about the axis as a function of the arclength 
s. In terms of the Euler angles, the rates of rotation are 


^3 


dO dcj) 

sm^p- -cos yjsmU-—, 

ds ds 

.dO . dcj) 

— cos yj- - smyjsmU-—, 

ds ds 

dp; dcj) 

—^-cost' — . 

ds ds 


To the lowest order in fl{s)^ the energy cost of ribbon deformations can be expressed as m 


(S43) 


E = Jy [^ 1 ^? + + cnl] . (S44) 

If ribbon is constructed from a 3-dimensional isotropic elastic material of Young’s modulus E and Poisson’s ration 
then the parameters Ai are [14] 


Ai = EWh^/12, A 2 = EW^h/12, C = 


(S45) 


where fis = E/2(l + z^) is the 3-dimensional shear modulus and u is the Poisson’s ratio. In terms of the two-dimensional 
graphene elastic parameters k, Y and u in the main text, we have 

Ai = A2 = YW^I12, C = 2 kW{1-u). (S46) 

In the limit of large Foppl-von Karman number YW‘^jn ^ 1, we find that A 2 ^ Ai,(7. As was shown in the 
main text, for ribbons whose width W is larger than the thermal length scale £th ^ hi/^/ksTY the internal thermal 
fluctuations of the ribbon renormalize bending and twisting rigidities to 


Ai ^ 




A 2 - 


C 


kW^+v 


(S47) 


where 77 ~ 0.82 and rju = 2 — 2r] ^ 0.36. As noted there for graphene membranes the thermal length at room 
temperature is of order the lattice constant ith ^ lA, so these renormalizations can be extremely large. 

In the presence of an external edge force E along the laboratory z-axis, the total free energy becomes 


E = Jy [^ 1^1 + ^ 2^1 + Crii] - Fz, 


(S48) 


where z = fds (es • z) is the ribbon end-to-end separation in the z direction. In the presence of thermal fluctuations, 
the expected value of ^ is 


(z) = keT 


dlnZ 
dF ’ 


(S49) 


where we introduced the partition function 


-E/keT 


Z = j V[&{s)]e 


(S50) 
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A. Schrodinger like equation 


By the usual transfer matrix/path integral arguments for statistical mechanics in one dimension, the partition 
function Z is closely related to the propagator for the probability distribution of ribbon frame orientation, where the 
unnormalized propagator is defined as 


G(©/,s/| 0 i)= f 

Je 


@(Sf) = ®f 


V[&{s)]e 


-E/ksT 


/®(0)=®i 

The function above propagates the initial distribution of Euler angles /9(©o,0) to 

/9(©,s) = y G(©,s|©o)/9(©o,0)d©o, 
where p(0, s) is unnormalized and the partition function is expressed as 

z = j d&p{&,L), 


(S51) 


(S52) 


(S53) 


with the Euler-angle measure fdS = dcj) sin OdO d'lp. In order to derive a differential equation for the 
propagator, we consider its evolution over a short ribbon segment Ss: 


+5s\&i) = jd® 


(S54) 


From the equation above we follow Ref. [25] to derive an imaginary time Schrodinger equation for the propagator 


)G = (5(s)5(©-©o). 


where H is the Hamiltonian defined as 


_kBT Jf Jl Ji\ F{es-z) 

2 [a, ^ A, ^aJ knT ■ 


(S55) 


(S56) 


Here the {Ji} are the angular momentum operators around the ribbon frame axes e^, which can be expressed in terms 
of derivatives with respect to Euler angles [24l |30|. The distribution of ribbon frame orientations obeys a similar 
differential equation 


p = 0, for s > 0. 

By expanding the distribution of initial ribbon frame orientation in eigen-distributions pa(0), where Hpa 
the ribbon frame orientation distribution and the partition function can be expressed as 

= y]Q!ae"^“®Pa(€)), 

a 

Z = [d&pa{&). 

a 



(S57) 


^aPa 5 


(S58) 


In the thermodynamic limit of very long ribbons (L ^ oo) the term with the smallest eigenvalue dominates in the 
partition function and the expected value for the end-to-end separation of the ribbon in z direction becomes 



(S59) 
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B. Analogy with the rotating top in quantum mechanics 


To proceed further (and to derive results valid for finite L as well as L ^ oo), we note that the differential equation 
(S57) looks like a quantum Schrodinger equation for a rotating top in gravitational field proportional to F, where 

Hence, we can borrow methods from quantum mechanics to find 


the coordinate s acts like imaginary time 
eigen-distributions pa and eigenvalues A^. For a rotating top it is convenient to expand eigen-distributions in the basis 
of Wigner D functions [24] with a well defined total angular momentum = J(J + 

and angular momentum projections along the ribbon tangent and the z axis JzDj^^iS) = 


oo J 


Pa{^) — ^ ^ ^ 


J=0 K=-J M=-J 


In order to evaluate the partition function Z in Eq. (S58), we need to evaluate integrals like 

oo J J 


[ d&Pai&) = E E E = 

J T _n _ T T\/r _ T 




a, 0,0 • 


(S60) 


(S61) 


J=0 K=-J M=-J 


Note that only those eigen-distributions Pq( 0), whi ch have non-zero component contribute to the partition 

function Z. Since the Hamiltonian H in Eq. ( |S56| ) does not mix Wigner D functions with different M quantum 
numbers [24|, we can restrict the search for eigen-distributions Pa(^) to the subspace with M = 0, where Wigner D 
matrices can be expressed in terms of the spherical harmonics 

e, 4>) = 0), (S62) 

where * denotes the complex conjugate. In order to avoid additional normalization factors, it is convenient to expand 
eigen-distributions in the basis of spherical harmonics 


Pa{tp,0,(t)) 


OO J 


E E 


(S63) 


Then the eigenvalues \a and corresponding eigen-distributions pa{ip,0,(j)) can be found from the matrix equation 


J2{{J'X'\H\J,K)-XSja'5K,K')C^ = 0, 


J,K 


where 


{f,X\H\J,K)= / sindde dcj) Y^'* {9,4>) H Yf {9,4>), 

Jo Jo 

{J{J+1)-K^) / 1 1 


{f,X\H\J,K) = 


knT . 


3j,j' 


tk,k'- 


— + — ) +dK,K'-^ 


As 




xJ{J + K){J + A - I)(J - A + I)(J - A + 2) / I 


F Sk,K‘ 


Ai A2 

xJ{J + K'){J + K' - I)(J - A' + I)(J - A' + 2) fj. _ 1_ 

Ai A2 


ksT ^(2J-|-I)(2J' + I) L 




(S64) 


(S65) 


We have solved the above matrix equation numerically to find the whole spectrum of eigenvalues and eigen- 
distributions pa('0,^,0). 

In order to evaluate the partition function E, we need to expand the initial ribbon orientation in terms of the 
eigen-distributions 


p27r 

P pi'll), 9,(j),S = 0)= dip p{lp, 6», S = 0) = W ttapaiP, 0, (p), 

JO 


(S66) 
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where P denotes projection to the M = 0 subspace and 

^277 


^277 



1*277 

PIT 


: / d'lp 

/ sin OdO 


Jo 

Jo 


p2tv 

PIT 

II 

-- / d'lp 

/ sin OdO 


Jo 

Jo 


oo J 


J=0 K=-J 


-iK^K 


The partition function is then 


(S67) 


rZTT PIT pZTT 

Z = y^aae~^^^ / d'lp sinOdO / dcj) pa{i^,0,(j)), 

a "'O Jo 


(S 68 ) 


Finally, the average ribbon end-to-end distance {z) is obtained by taking derivative of this partition function Z with 
respect to force, see Eq. (S49). 

As was mentioned in the main text, this same formalism can be used to study both the pulling and bending of 
ribbons. For pulling we orient the ribbon along the z axis with the initial ribbon orientation = {0,0,0}, which 
results in 


=Yf*iO = 0,^ = 0)= 5k fi (S69) 

For bending around axis ei we orient the ribbon along the x axis with the initial ribbon orientation 0^ = {7r/2,7r/2,0}, 
which results in 




= F/* (0 = |, 0 = I) = 


2l^l(-f)^ 

- COS 

w+im 


27r 

2 

V 


'(2J +l)(J-|7f|)! r[(J+|7f| + l)/2] 
(J+|i^|)! r[(J-|i^|+ 2 )/ 2 ]' 

(S70) 

For bending around axis e 2 , which is harder because it involves the ribbon stretching, we orient the ribbon with the 
initial orientation 0^ = { 0 , 7 r/ 2 , 0 }, which results in 


= F/* 


2 l^l 

pV + 1 ^ 1)1 

h 

-cos 

27r 

2 

r 


/(2J + l)(J-|7f|)! r[(J+|7f| + l)/2] 
(J+|i^|)! r[(J-|i^|+ 2 )/ 2 ]- 


(S71) 


C. Large force limit 

For large pulling forces, we have to take into account both the ribbon stretching and the deformation energies that 
appear in 


E = J^[Ajnl + A2nl + cnl + FiduL] - IdsF{z- 63 ) [1 + u,,], (S72) 

where Uzz corresponds to the stretching strain along the ribbon backbone, and the one dimensional Young’s modulus 
is Yid = Yr{W)W. Here, Yr{W) is the renormalized 2 -dimensional Young’s modulus evaluated at the scale of the 
ribbon width. For large pulling forces the ribbon is nearly straight and the tangent es can be approximated as 


eS = + tyf + 



(S73) 


where ^ 1. To quadratic order in tx and ty , the free energy becomes 



T 



dt 

ds 


- tyfts 


C o Yid 2 / 2 2 \ 

+ -^^3 + ~^‘^zz + -^ (4 + ty) 


Fu, 


(S74) 
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After integrating out the Qs and Uzz the effective free energy becomes 


ds \ , f dt 


-r\ A 


+ ^2 ^ + F (tl + tl) - 


2 , [Ait^{dty/ds) - A2ty{dt^/ds)Y 


For C > 0 the last term is 4-th order in and ty and can thus be neglected for large forces. Upon rewriting the 
effective free energy in Fourier space 


^ E [(^ + + {F + A^q^)\ty{q 




L{F + A2q^) 




Using the results above we can find the ribbon extension 


ds (z • 63 ) [1 +' 


ds 1 + - 


j) = 1 + K.)-+ 


L/ Yid 4a/U \a/Ai a/A2 / 


We find that for large forces the ribbon extension is indepe ndent of the twisting rigidity C. 



